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The propagation and mutual transformation. of acoustic and 
thermal waves is investigated in media with a high radiative press- 
ure, The investigation is carried out by means of the equations of 
hydrodynamics for m,atter and of the radiative transfer equations in 
a moving medium (in the first order with respect to v/c) in the 
Eddington approximation. A number of model problems which explain 
the physics of the phenomenon are solved in the 'first part of the 
.paper; Waves in a homogeneous medium with an abrupt jump in opacity 
and in a medium of variable opacity. In the second part the outflow 
of waves is studied in the 'atmospheres of astrophysical objects; 
the accretion disks around black holes and very massive stars. The 
characteristics of the luminosity variability 'of Cygnus X-1 and of 
the nuclei of galaxies are examined in- these models. The presence 
of convection and turbulence generates acoustic. waves whose spectrum 
upon outflow inco transparent layers is dete,rmined by the conditions 
of transmission and of damping. Variability of t^ie radiation is 
connected with fluctuations of the photospheric' temperature and that 
of the corona because of variable heating. The characteristic times 
•of variability agree well with observations for; all the objects, but 
it is difficult to obtain a sufficient amplitude for the brightness 
fluctuations for very massive stars. . ■ ■ ' 

I. THE BASIC EQUATIONS AND THE CASE OP A HOMOGENEOUS MEDIUTVI 

INTRODUCTION ’ ' ' 

• f 

In astrophysical objects conditions are often encountered where 
the radiation pressure significantly ex.ceeds the ’gas pressure, 
although the material mass is determined by the, gas, • A typical 
example are m.assive stars, which are well described by an Emden poly^ 
trope with n=3, and where the ratio 3 of the gas 'pressure to the 
radiative pressure equals (see [1]) 
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Here y is the number of nucleons for one particle; for hydrogen y=0.5. 

f 

At M>72Mgj the radiation pressure is already greater than the gas 
pressurej and for very massive stars with M>10^M^, which are used as . 
models for a quasar or for the nucleus of a galaxy (see [2]), g<<l. 

For these objects and also in the theoretical investigation of x-ray 
sources, a model of disk accretion onto a black kole is often con- 
sidered (see [ 3 ] through [6]). For a fairly large luminosity L, 

3<<1 occurs for the central region of the disk. 

Convective instability is a characteristic property of regions . ■ 
with a predominance of radiative pressure ■(3<<1)» For very massive 
stars the presence of convection follows from the .condition of radia- • 
tive equilibrium (see [2] and albo [7l)> and convective instability . . 
of regions with 3<<1 for disk accretion were shown in [B] (also see 

[9]). . ... 


Convection in the subphotospheric layers of a star or of a disk • 
will i-inavoidably lead to the appearance of a mechanical, undulating 
energy flux which disappears at a small optical thickness x<<l, being' 
transformed into heat. This leads to the formation of a hot gaseous ‘/6 
corona that is similar to the solar corona, whose temperature T^ is 
much greater than the photospherlc temperature ' T^^. One succeeds in 
'explaining certain features of the Cygnus X-1 radiation source with . 
the presence of a corona in an accretion disk around a black hole 
(See [8], [10] and [11]; moreover, the structure of the corona can be; 
determined to a significant degree by the ma^etic . field (see [12]). 

■ It is also possible that the examination of the -coronas around very 
massive stars will turn out to be necessary to. explain the properties- ‘ 
of certain quasars or active nuclei. ’■ 

The flux of mechanical energy that is generated in the convective- 
zone is transferred outward, by sound or Alfven and magnetosonic 
waves into layers with r<<l through the outer 'radiative regions. If 
the radiation pressuriy is negligibly small in 'the convective zone and 
in the photosphere, then the mechanical energy flux (with no magnetic . 

2 ' . 


field) is where v is the velocity of the material in the wave and v^ 
is the speed of sound that occurs practically without damping in a 
region where t< 1. Transformation of the wave into a shock wave and 
dissipation of the magnetic field occur there, and the mechanical' 
and magnetic energies go into heating the corona. If 3<<I> then only 
a small fraction of the flux P that is generated in an. optically 

ITl . . 

thick region goes to beat the corona (see [8]). The great attentua- 
tion of the mechanical energy flux upon its outflow into a trans- 
parent region under the conditions 3<<i is connected with two circum-' 
stances. First, upon going over from a region where t>> 1 to one 
where x<<l, the speed of sound which, inside this zone is related 
to the radiative pressure and outside of it, at t<< 1, to the gas 
pressure, decreases: v^ (x<<l)/v (x>>l) - 3“/ 'The second reason 

for the decrease of P^^ is the strong damping o,f the sound waves that . /7 
is connected with radiative friction and with .radiative heat con- 
ductivity. . , . ■ 

A crude estimate for the total attenuation ' of the energy flux • 

by 3”^ times was used in [8]. The propagation, of sound and thermal . ■ 
waves in a medium with 3<<1 and the transfer of mechanical energy 
into the transparent layers' above the photosphere- are investigated 
in the present paper. .The entire inve'stigation is carried out on the’ . 
basis of the equations of hydrodynamics for matter and of the equa- _ • 

tion of transfer for radiation, which is used ;in the Eddington approx- 
imation which allows one to describe the transparent and opaque . 

regions in a single form. A flat geometry is ’Considered and the 
effects of sphericity are neglected. -Terms of the order of ~ v/c 
are taken into account in the equation of transfer as is the case 
also in [13] and [l4]. ^ 

In the first part of the paper, a simple , derivation of the • ‘ 
equation of transfer in a moving medium is given' 'in Paragraph 2 arid 
the basic equations are shown there-^ ’In the following paragraphs 
the propagation of waves is Investigated without a magnetic field 
.under the conditions 3<<1: of a wave in a homogeneous medium and ; 

in a medium with an abrupt jump of opacity. In the second part the 
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outflow of a wave from a plane, static -atmosphere is invesgitated. 
Numerical estimates are. made for the flux of 'mechanical energy in a 
model of disk accretion and in a very massive star. The estimate 
of a damping of 3”^ times that has been used in [8] turned out to be •. 
fairly good. The condition for the outflow of a wave into a trans- 
parent region singles out a characteristic frequency whl '’-h can be 

.* * * 

connected with the observed frequencies of the fluctuations and 
brightness variabilities in the x-ray sources Cygnus X-1 and 
Gircinus X-1, and also in some nuclei of galaxies -and in quasars. 

In conclusion, numerical estimates are carried -out that are connected 
with the existence of the characteristic frequency, and the observa- 
tional consequences of the model are given, ' 

I , Set-up of the problem and the basic equations . _ ^ 

a) The e.Quayon_of ;tr_ans_fer_in a. moving_m^dium 


Upon’ the outflow of a wave into a transparent region, it is 
necessary to use an equation of transfer with allowance for motion 
of the matter. In regions with x>>l it is sufficient to limit one- 
self to the equations of radiative thermal conductivity; however, 
for T<<1, the flux and 'the density of the radiative energy are not 
connected unambiguously and it is necessary to consider at least the’ 
Eddington two-moment approximation. The equation, of transfer in 
•moving matter with allowance for ~v/c was derived by a fairly complex 
"geometric" method in [13] and [l4]. Here we shall show a simpler 
method for deriving these equations that is mentioned in [Ijj but. has 
not been carried out there. In an arbitrary inertial frame of refer- 
ence and in the framework of the special theory, of relativity, the 
equation of transfer' has the form (see [15]) 


c Dt. 




( 2 ) 


Here lu^ is the spectral intensity of the radiation, is the fre- 
quency, Jq and k^ are the coefficients of emission and of absorption, 
which include scattering, and n^ is the direction of propagation of , 
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the radiation. All the quantities with an ir>dex "o" are determined 

I 

in the inertial frame of reference. The coefficients of emission 
and absorption are most simply express, ed by the properties of the 
matter i.' the system with respect to which the 'matter is at rest. 

Let us call this system a Lagrangian system. Let us limit ourselves 
to a plane, one-dimensional case where all' the 'quantitiqs depend 

Then 


only on the one spatial variable z^. 


-V, = 


. ^3 

■) 






co*^ t n, V. ) , 


( 3 ) 


If the matter is moving with a velocity v in the z-diroction, then 
the frequency o and the cosine y of the angle in the Lagrangian 
system are connected with the quantities in the inertial system by 
relations (see [1^]) that are written down with an accuracy ~v/c: 


n) (i- V (<-■» 

|a=|a. -W c u. r ^ ^ (V v/c. 


(^) 


The distribution function n(p, a, t) of the .photons in phase 


space is related to the intensity by the relation (see [16]) 


\\ rr --- ~~ 


( 5 ) 


Prom the invariance of the .distribution function 'n with respect to 
Lorentz transformations, invariance follows for the quantity 


L.K = uh\ 


( 6 ) 


Using Equation (6), one can write down the left-h'and part of Equa- 
tion (2) in an invariant form and then the invariance of the right- 

% 

hand part of Equation (2) will be reduced to the relations 




( 7 ) 


•fhe equation dt^/u^=dt/u was taken into account here and Relations 


(6) and (7) were used in [13] and [l4.]. 
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Let us write down Equation (2) in, the Cm,’t) Lagrar>{?ian system 
where, with the necessary accuracy with respedt to v/c . 




JpCeOdl’/. 

5, 


( 8 ) 
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Taking the transformations and ip:variants (^')» (8) and (7) 
into account, Equation (2) with allowance for Equation (3) to an 
accuracy of v/c, is written down in the form 


A '^vl 

c <'d-t 


/ c, X ■T.it 'ijv 


(9) 


'1L> ^ _ y. 
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One can write down the equation of continuity -in Lagrangian coordi- 
nates in the form dz/dm = 1/p or, after differentiation, 

^ ; (10) 
After substituting Equation (10) into Equation (9)» we shall finally' 
obtain an equation which also follows from [1^] for r-»-“ ; 


c H'’ ^ I 




(11) 
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If, in addition to pure absorption, scattering occurs, then one 
can, for a condition of local thermodynamic equilibriam (LTE) , write' 
down (see [17]) ' . ’ 

- oc,e>AT) ^ cs, ji,ciQ/^:ifr, . ■’ 




(12) 


6-^ = G.65-i0"^^ CJM». 

Here au and a are the coefficients .of • absorption and scattering for 
electrons, B (T) is the Planck function, n is the concentration of 
electrons, and is the Thomson cross-section. ’.The scattering is 
considered to be coherent and isotropic. . These approximatons are 
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sufficiGHt ly good for the conditions being considered. 

b) Equations^ for the JBoments^ in^the J 2 ddin£ton_A 2 proxim^ Zll 

A system of angular moments that is obtained after multiplying 
Equation (11) by and integration with respect to y from -1 to 
+1 is often used for an approximate solution of the problem instead 
of Equation (11). The first two moment' equations have the form 

C J r^tr^f ''?V ■ ■ 


4 i-.:|C(3K.- J.v c-.a^Cr) 




cot J Ovn . cp^t'OV .<p '3’^ • ' 


— — ^ ^ V ' 

The 'following moments are introduced here: 


' 5 i i j 1^' F * • 

-V ' • . 


(14) 


If one is not interested in effects connected with lines, then 
one can average Equation (13) with respect to frequency. We obtain 


i;®? 4 p9.h! 


=. c^p boO 


►o ^ 

- j O^v Jv *=iw. 


(15) 




(16) 


Here a is the absorption coefficient averaged with respect to the /3^ 
Planck^ function B(T) = acTV4iT, a- 7 . 56 • lO”^^ m .cgs’ units'. The average 
absorption coefficients with the weights and H^ are usually' expressed 
4 by Cp and the Rosseland mean, .respectively (see [14]) 
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i ^\> Jv B- C<p ’7 , J H^®iv sf 

J^JJvJv, H-jVl^<-|v>, 


(17) 


In later calculations we shall consider as independent of u (the 


grey approximation) and Oip=ap^*= a. Let us ‘notice that, upon averag- 
ing with respect to frequency, the third moment falls out .of 
Equations (15) and (l6) for the moments. 


We shall solve the problem in the Eddington approximation, that 
is, let us set 


K = 4 J 


(IB) 


This approximation is fulfilled with high, accuracy for a large 
optical thickness t, and for small t values the error is small; we 
have K=0.'i)lJ in a plane, grey atmosphere at t=0 (see [17] and [18]). 
Finally, in the Eddington approximation, the e^quatlons for the 
moments have the form: • ' ’ . 




c 'd-t 


c) Hydrod;^namic_eq_uatrons_f£rjnat;t^ ’ •. . ; 


- p 

Z'SX 0)0'^ »vv CjO' 


(19) 

(20) 


In a gravity field with a gravitational acceleration g=const>0 
the equation of motion has the form: . . ' 


Q)m 


roi> 


3 . 


■Qi 


( 21 ) 


Here the gas pressure is P --RpT/y. R=8.3*10 in cgs units, and p is 

S 


the number of baryons for one particle.. The last, term in Equation 
(21) represents the force acting on the matter' from the direction of. 
the radiation flux. The equation which takes into account the ex- 
change of energy between matter and radiation, has the form: 


8 
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( 22 ) 



where Eg=RT/y(Y“l) and y is the ratio of specific heat for the gas. 
One must add the equation of continuity in Lagrangian coordinates 
Equation (10) to Equations (21) and (22). We shall show that the 
complete system of Equation (10) and Equations (19) through (22) 
leads to the correct form of the law of conservation of energy with 

an accuracy to terms '^(v/c). The density E„, ergs/gm, pressure P„, 

2 2 ^ ^ 
dynes/cm , and energy flux ergs/cm »sec, of the radiation are 

expressed by the moments J and H; , ■ ‘ 


E, = _^c Jiv civ -.p jlv c'Hv'.--. , (23) 

R-ipE.- J. ' ' ■ ;• 

E jl„ j. JQ jv = Js) - ■'iJrH . 


Multiplying Equation (21) by v, Equation . (19) by Hir/p and adding 
this equation to Equation (22), using ’Sv/Sm for' the expression 9p/9t 
in Equation (10), we obtain 







dc qim 


J 








(2i|) 


Let us express the last term in Equation (2^}’) -by means of Equation 

2 

(20), Neglecting terms of higher 'order than (v/pc )3H/3t and .Ziii 

(vH/c )3v/3m, we obtain, taking Equation (23) into account, an equa- 
tion of the conservation of energy in’ the., form: , 




(25) 


Here on the left stand the thermal, radiation',' potential and kinetic_ 
energy for a unit of material mass under 'the 'derivation sign, the 
first term on the right expresses the' work of. the pressure forces, 
and the second term is the divergence of the radiation flux. 
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# • 

I 

Xt is more convenient to introduce, the opacities 

a. « Cc^■^ ^)/p , cx/p . • (26) 

in Equations (19) through (22) instead of the coefficients a ar,id a. 


Under the conditions considered bplow, the cosificient of brak- , 

I, 

ing absorption a is much less than the coefficieht of scattering for . 
Electrons (see Cl8]): 


, ■> "T " * ' 

^ 6 . 6 - 10 '*® rv 


(27) 


On the other hand, at sufficiently high temperatures the energy • 
of a photon changes due to Compton processes (fqr an invariant num- 
ber of photons) , which will lead to an exchange of energy between the 
gas and radiation and must be taken into account in. Equations (19) 
and (22). In view of the great complexity of a rigorous allowance 
for non-coherent Compton scattering in a non-homogeneous medium, we 
shall limit ourselves to two limiting cases. In the first case we 
shall completely neglect ab.sorptlon (pure scattering) ^ ; 


(I) 


- Gyp . ae, . 0.2 c i •« Xii) <:.uV'‘-, 
ce, - 0, 


P Vl5 

cm*/gm, ( 28 ) 


' Xpj is the concentration of hydrogen by weight. 

In the second case lob us consider Compton scattering which 
acts as a source of absorption: ■ ’ • 

(II) . oe, . ■ 

p ■ (29) 

’.In the examples considered below the fraction 'of 'the mechanical 
energy flowing out from the photosphere is practically identical in . 
Cases (I) and (II). Actually, the case of Compton scattering corres- 
ponds to an intermediate case 0 < < a:-, , ■ . 


, . ilirfr -- 
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2 , Thfi llnoarlgQd eauatioAs for n mll disturbances at 

Let us introduce the Lasransians of the small disturbances 

0«ST, j «. ^20) 

with respect to the equilibrium solution of Equations (10) and (19) 
through (22) : 

P Crv») , •) ^ ) H *■ con si , 

'iCvn*) ^ B" acT^/^/a». ( 3 ;L) 

Let us also set all the quantities of Equations ;( 30) a exp(-iwt)i 

- - - 4 ^ \ j. < \ 4- Q VH • 


- - h 

cin ^ ' 9' ^ . 

(32) 

;o* ae,Cl- ^ 

c-J.^-'ann 30r ‘^0 A 

(33) 


(3^0 

» ^ t r* ' ^ X2L^* f-v . 

(35) 

-l-i • 

(36) 


Eliminating p from Equations (33) through (36)_ by means of Equation 
032), we shall obtain • ^ • ' ■ _ ’ 


ihii - dii ^ qJ- ;o d|j ~ ^ o_c-.-sa\ 

CpJ ciin ^ 30T WWv af* 


-<» se,acT^ ' 


(37) ■ 

(38) . 

(39) : 

(itO) 
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In Equation (28), Case (I) , setting we have from Equa- 
tions and (37) for Case (I): ' , 

O'* ) ‘ (^la) 


H . I . A r«^ * ' 




33T 


Mm 


Cjf» < 

♦ ... 

In Equation (29), Case (II), let us immediately allow for the 


smallness of the quantity $®Pg/Pj^«l, and also' let us consider the 
'case where the inequity 

/w.pc << i . 


(il3) 


is fulfilled. This inequality is fulfilled at oscillation fre- 
quencies that are not too high. Taking inequality (^3) into account 
one can neglect the second term in the denominator of Equation (kO) 
and obtain, for Case (II); 

ft. 

'i*. <vc-r“ ■ (i(ib) 

Taking Equation (4lto) into acoour.it,, the quantity in the right-hand 
part of Equation (37) equals ^ 

tc.Cj-«.cT''eA) = - ■ (b) 

By virtue of the inequality 8<<1, this term is much smaller than the 
last term in the left-hand part of Equation (37) and one can neglect 
it. Thus, Equation (k2) is valid both in Case' (T) and also in Case 

, ■ f 

(II). Let us now eliminate 6 from Equation (39) • .Using Equation 
(^lla), in Case (I) we obtain from Equation (39.)' • 


/ ' ■ (nila) 

In case (II), using Equation (4lb), we obtain ,from Equation (39), 
with allowance for inequality (43)s 

■ ..Xu - -oO f 

c . . (44b) 

Equations (38), (42) and (44) will be used in .lafer calculations in 
regard to the variables y, h and j. 


3. Waves in a homogeneous medium ^ 


, *This case was also considered in [29]. 
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In this case H»0 all the equilibrium va;Lues in Equations (31) 
.are constant, among them p»const, and instead of Inequality (43) 
one can use the stronger inequality 

to /c >1 Qjo ( 45 ) 

This inequality indicates the smallness of the free path .of a photon 
d s in comparison with the distance divided by 2 it which a 

photon pisses through freely in the period of the vibrations. One 
ean look for the solution of Equations ( 38 ), (-42) 'and (44) in' a 
homogeneous medium in the form exp (ikz). . Als’‘o taking Condition (45) 
into account to simplify Equation (38), we- shall obtain 


(46) 



k . - 


= 0. 





!?) 







i 0 

cp 


sJ." 

r ^ ■ 


(I) 

/f jt *5Co 

c 

1 

(o' - 

f) 

^ P 



(II)' 

c 

' k - 

d 2T. 
'^P 

-h 

■* ‘ 

p 


Equating the determinant of the system of 'Equations (46) to zero, 
introducing the notations ' ^ 

,..v ii Pr r-i _ . n ■ J._ 

"c ’ , 5) " ^ f ’ ^ ' (/I7) 

and allowing for the fact that 3<<1, ye obtain- the dispersion equa- 
tion in the form: 




(M) 

For Case (I) the adiabatic speed of sound Vg from Equation (47) 

enters into Equation (48), and for Case (II) the Isothermal speed 

2 ’ ' 
of sound Vm=P /p enters instead of v„, that .is', 'Equations (47) and 

K w , 

(48) are applicable for Case (II) if one sets y-1. 

In this paper only forced vibrations (accprdlng. to the termino- 
logy of [ 19 ]) are considered, where co -is considered real and given, 

p 

and k is found from Equation (48). ' Ife obtain 
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• ' ' 

( 49 ) 


A ' ^ i. . 


Taking the 
we obtain 

2 2 ■ ■ • 
smallness of v /v '^3<<1 into account in Equation 
B ^ 

( 49 ), 

< 

A* ^ 4 • ato • 

(50) 


pa to'*- ■ to=* ■ . wc t \ 

(51) 


The solution of Equation (50) describes the prqpagation of waves 
that are associated with the gas pressure. If absprption is absent 
(Case I), then the waves are adiabatic, and in ’the' case of strong 
coupling of gas and radiation (Case II). the speed of sound becomes 
isothermal. The solution of Equation (51) desdribes adiabatic waves 
that are associated -with radiation pressure. The frequency 

(52) 
(53a) 

(53b) 

(54a) 
(54b) 

It follows from Equation (53) that the low-frequency distur- 
bances of the radiation pressure are propagated with a velocity 
v^ (Equation (53b)) almost without damping, and the propagation of 
low-frequency disturbances in the gas assumes th'e character of 

14 ■ 


is critical for both types of waves. Eor w<<y^, we have 


For o)>w^, we obtain 

t § (i - • 
k ^ -L ^ 4*0 A “ ' — Y 










"diffusion" waves (Equation (53a)) with a coefficient of diffusion 

2 ’ ‘ 

. "’hese waves describe the leveling of the inhomo- 

■geneity in the gas ("entropy" disturbances'). On the contrary, the 
high-frequency disturbances in the gas are propagated as sound waves /20 
(Equation (5^a)) and in the radiation they assume, the character of 
thermal waves (Equation (5^b)) with a coefficient of temperature con- 
ductivity “1/3 q 1 (compare with [20]). • Let . us 'remember that, for ' 

. the waves in a gas are isothermal. 

A medium with an abrupt ,1ump of opacity '. 


Let us consider a medium in Which p, T and are homogeneous 
and the opacity is . It undergoes an abrup.t jump in the plane 
s=0. The quantities 1 and to-j^ undergo the same kind of abrupt jump 
in this plane. Let '^e. C le. o’) » at. A wave which is propagated 

in the direction of increasing s can change the nature of its propa- 
gation at B=0 if A • /> ' 


« O C< 

a. ' - 




b 


pw.a<p') ' (55) 

Short-wavelength disturbances in a homogeneous medium w>>w^ rapidly •. 
damp out even in the case where their propagation has the nature of 
the wave of Equation (5^a). The damping length equals l^^-2v^/u)2= 
2/3cv 1/v^p, which at small. 3 can become, on the -order of the length- 

o 3 X . 

of the free path 1 of the photons. On the other hand, the damping 

P • 

of long waves is very weak; .from Equation (53b) we have 
=6v^/clw->-“ for w->0. Thus, it makes sense to investigate the transly 
tion of only the fairly long waves through the abrupt jump at to=0. 

The case of Inequality (55) that is considere'd below is the most 
interesting one from the point of vle,w of later applications. 


The wave of Equation (53b), upon; passing through an abrupt jump- 

of opacity, is split into four waves:’ the two_ reflected waves of 

Equation (53) and the two transmission waves which, for s>0, a)>t 0 g 

already belong to waves of the Equation -(5^) .type. ’ It is obvious that, 

at the point of the abrupt jump of opacity, the y displacements, dis-. 

turbances of the gas pressure P =P (p/p -te/T) and of the first moment /21 

s s • * 

'■ ■ ■• 15 ■ ■ 











J which characterizes the radiation pressure of Equation (23) remain 
continuous. The continuity of the second moment -h also follows from 
the condition of the absence of an energy sink or source at the 
place of the abrupt jump. If is the amplitude of the incident 
wave, and A^ are the amplitudes of the reflected waves of the 
Equation (53) type, and and Bg are the amplitudes of the trans- 
mission waves of the Equation (5^) tyr then one can write down the 
continuity condition for the physical syjr.n titles; at the abrupt jump 
at a=0 in the form 


^ ^ ^ 1 il ■» • 


(56) 


To obtain the amplitudes of A^,^ Ag, B^ and Bg ,as functions of A^,' 
one must express the amplitude',? for different q values by some 
defined q^, using Equations (46), (32) and (4l), and solve the four 
linear non-homogeneous equations which have been obtained. If one 
adopts Qg=j j then one can write down the remaining continuous quan- 
tities, using Conditions (45) and 3<<1, in the' form: 


a _ _ : ; a ^ A; _ h ■ 

y sac.yoy "it. > 

(I) S -- . (3 ^ ^ , it) P? , (^' * ■; ^ 

p r.cl ^ pg e-cT 


(57) 


The solution of the system of Equation (56), taking Equation (53), 
(54) and (57) into account,’ has^ the form for q=j , • 



(58) 
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For q=y we have 


A.;,- -zyi^c-O ['- t>-0lro, - C‘-0Nfl 1 A,j‘ 

Av^ ■ A.J, 


B.j,= 

b.a= A.,'. 


(59) 


The relations ( 58 ) and (59) ai*© valid for both Cases (I) and 

(II) if one assumes v to be equal to the adiabatic and the isother- • 

S 

mal speeds of sound respectively. The term.s and v^/v^, 

were neglected in comparison with unity in all quantities. It follows 
from Equations ( 58 ) and (59) that, upon passage through an abrupt 
jump of the type of Inequality (55), an almost -complete reflection 
occurs for a radiation wave of the Equation (53b)'_type. The ampli- ; 
tudes of the remaining j -waves are significantly smaller than ^ 2 ^'* 
the shift of the discontinuity point generates a transmission gas 
wave of almost the same amplitude 


I Bijj 1 J. 


( 60 ) /23 


Thus, upon passage through the abrupt jump of opacity, a wave trans-.. 
formation occurs such that the wave which 'possesses the minimum 
damping "survives". Relation (60) shows that y, the amplitude of 
the wave v;ith minimum damping is almost continuous at the abrupt jump. 


We also considered the question of a wave transformation for the 
case where the opacity (and the paths of the photons; t ) change con- • 
tinuously. Here the system of Equations (38), (39) and (^2) is 
reduced to :one equation for the displacement of, y'(a dashed line 
indicates a derivative with respect to a) i ‘ ' ■ 
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It was assumed that for a ^-6 the path g is oonstant , for 

•2:’>§ -^-^5 > 6<5 v ^jQthe path increased- monotonically from' 

to , and, moreover, the relation C.(?) was selected in the 
form of a fifth degree polynomial in such a manner that • 

t' c< (^^ ' (62) 

Then , and are continuous for all a. -The parameter the 

half -width of the transition region, varied in different variants 
of the calculation from values to , where 

■ ( 53 ) 

is' the wavelength in the gas (without allowance for damping). 


Equation (61) was solved numerically by an' implicit fifth order 
Adams type method (see [ 21 ])’. In order to obtain a solution for 
Equation ’(6l) which corresponds physically to a radiation wave going 
from the region 2 ^ & > we started in the following manner. First,, 

a "radiation" wave going in the direction of increasing a was given 
at ‘i='i t . From this the initial conditions for y'. (that is, y, y’ , ’ 

y*' and y''*) were obtained for and an integration "backward" 

of Equation (6I) was carried out from 2 = to , From the /24 

y, y', y’ ' and y’ ' ' values at ■£---£ the amplitudes of four waves 
(radiation and gas, incident and reflected) were found. Next the same 
procedure was repeated where the initial conditions for corres- 

ponded to a purely "gas" wave. A linear combination of these two 
solutions in which the amplitude of the incident gas wave equals zero 
for -2= -6 is' also the solution sought (a radiation wave is incident; 
waves of both types are reflected and transmitted). To check the 
accuracy, a repeated calculation was carried out from 2--C to IS 
with the initial condition found for . ‘ The 'accuracy is usually 

better than 10 ~^ and only for is an error of =^10“^ accumulated. 

The results of calculating amplitudes are shown in Table I for 
values of the parameters . ' . 
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fe,/fA I0”2, « Ib*^, 

Uj/u?^ « 10“^, co/tOjj^clO, 

C/>^ »o.oi, o.i, I H 10. 


(6il) 


It is evident from Table I that, for small val.ues of . the 

results agree well with Equation (59) that have been obtained for 
an abrupt jump of opacity (upon comparison it Is necessary to men- 
,tion that the quantities discarded in Equation (59) amount to 
in this example). Even in the case of Gs'X„one can use' the relations' 
for the modulus of the amplitude for an abrupf-jump by an order of 
magnitude (this case is most important for us 'for later applications)'. 
And only in the case of an extended transition zone 6 « 10 are the 
amplitudes of the transmission gas waves smaller than would follow 
from the simple principle for the transmission of a wave with a 
minimum damping. Here the point is t|iat the gas waves are formed 
in a region that is larger than their- wave length, they are coherent • 
and, by interfering, they extinguish each other (Ya. B. Zel'dovlch 
turned our attention to this fact). 


II. ASTROPHYSICAL APPLICATIONS 
5. A plane equilibrium at'mosphere 
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. Let us now go over to considering the propagation of waves in 
astrophysical atmospheric objects with 3<<1.- We shall use the 
approximation of a plane atmosphere in a constant field of gravity 
which, in the majority of practical cases, allows one to obtain 
sufficient accuracy. Equations (10) and (19) 'through (22) for a ■ 
static atmosphere assume the form , , . 


^ _ A_ 

5> ' 




H, 



At; ar . 


(65) 


The solution of the system of Equation (-65) for a known steady flux 
ha. the form 
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r - %-■ -Si. ' '"•’'•f"-’"'). 


( 66 ) 


For H<<H , the solution, Equation (66)^ is redu ed. to the well- 

C 

known solution for a plane atmosphere (see [17]). For the solution, 
Equation (66), the characteristic frequency .is,‘ from Equation 
(52), a function of the optical thickness t: 




( 67 ) 


If a wave is generated at t^>> 1 and its frequency m<w^(T^), 
then, at a)>to^(0), there exists such a at which ,w=w^(t) , and, 
at T<T^, the nature of the wave's propagation changes. At t>^ 1 
z<0, the solution. Equation (66), has the form • • 

t'I 

- 

Utv 




( 68 ) 


For T<<1, a>0, from Equation (66) we have: 


JT-P). o t'*- At'’- 


( 69 ) 
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7i = -^.Cvx\x' for 

With such a choice of the constants of integration the quantities T, 
p, and 5 are continuous for a transition from ;Equat ion (68) to 
Equation (69) at the point x=2/3. 


6 . The propagation of waves In a plane atmosphere 


The waves in a plane atmosphere are described by Equations (38), 
(j^2) and (^M. Going over to the variable s instead of m, let us 
write them down in the form 


vio : _ ^ + Sik \io 
c i 


ciz 

»< 

\iO ' 

30r 


l-l =-0 

se.c. > 


(70) 
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^ a) The ca^e__of a large optijcal thicknes£ 

The equilibrium solution is defined by the system of Equation 
(68). One can obtain a solution of the system . of Equation (70) in 
a quasi-classical approximation • . ' 


j - ^('^) (-V . ; 


(71) 

Leaving the main terms of the expansion, we obtain, by substituting 
Equation (71) in Equation (70), ' . ' 


P 0 ^' cj.pfv " ■* = 0 




30C, ■ 
S 


(72) 


(I) 


Let us use an inequality of the type of Inequaliti’es (i|5) 

t£.p„ 1 • ■ . (73) 

with from Equation (68). The system of Equation (72), which 
differs from the system of Equation (i|6) describing the propagation 
of waves in a homogeneous medium, only by the presence of the last 
term in the first equation. If one uses' the inequalities 

^ and /fecT c< 1 ^ 

which are assiamed as being always fulfilled, then -the dispersion . 
equation is reduced to Equation (48) if the quantities v , v^ and 
£= are considered to be variables defined. by the system of 
Equation (68). ‘ • 

21 


727 


maCASs-nmem^ 


As was noted in paragraph 5i the only t^pe. of waves that propa- ; 
gate with low damping in an optically thick medium are low-frequency 
(w«o)^) waves that are defined by Equation (53b)./ The damping of 
this wave is given by its imaginary part which, in the case of 
Equation (71), has the form ■ /2’8 


ie.garar/H ‘ ' 


(75) 


The amplitude of the displacement A, which equals A^ at t=t^, has 


the form 




(76) 


To calculate the factor f(t) in front of the exponential, one can • 
use the following expansion terms after substituting Equation (71) 
into Equation (70). However, this requires fairly cumbersome cal- 
culations. Therefore, let us find f(r) from simple physical consider'.- 
ations. Along with the wave damping, which is contained in the 
exponent's index, the amplitude of the wave must increase during its 
propagation outward through decreasing density. This intensification 
must also be exactly contained in- the factor f(T); If there were 
no damping, then the energy flux P transferred. by -the wave would 

be constant. ^ o V ^ _i . / it, x ‘ 

^ COhst . • . 

Taking p from Equation (68) and k from. Equation ’ (53b) into account, ; 

. F~ fcv)-x-c‘. ; • : 

Prom Equations (78) and (76) we obtain the change of the wave-'-s 
amplitude at 0 )<<w^ in the form 


mU'X <i'T eO , 




The value of Ti up to which Equation (79) is ^^alidv equals, with 
allowance for Equations (6?) and (68); 
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CCV,yOi» ^ ’’ 


( 80 ) 


V, 


vU 




In the majority of cases we shall consider waves for which T;>1, 
then, from Equation (80) it'follows that 

, ( 81 ) 

Which can be fulfilled simultaneously with Inequalities (73) and (74) 
If Ti>l then, at T<Tithe wave changes its character and "begins to be 
determined by Relations (54). Although a gas wave (Equation- (54a) ) 
is dampled slightly after one vibration period; for a given length 
Az it is, for t=ti, damped much more strongly than a thermal wave 
(see Equation (54b)): 


;»> i 


( 82 ) 


The amplitude of a thermal wave (Equation (54b)) in the quasi-classl-_ 
cal approximation of Equation (71) has, using Condition (77) to ■ 
determine f(x) and the solution of Equation (8), the form: 




(83) 


ma;*" 'c c 'v, . 


The damping for gas and thermal waves in Equation' (54) becomes ident- 
ical for . . ./ -1 ^ 

Im 4.. (84) /30 

If T2<1 in Equation (8^1), then the solution Equatdon (83) is extended 
to T=1 after which one must use a solution for an optically thin 
region. Otherwise when, - * . : , ' 


\ 
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c \r ^ • 
3 




(85) 


there is a region Kt<T 2 where a gas wave (Equation (5^a)) has less 
damping. Using Equations (5^a), (71), (77) and .(68), we obtain 




( 86 ) 




b. The c_as_e__of ^mall optical bhicknes_s 


In the limit of small optical thickness T<<l'the solution, 
Equation (69), of Equatlbh ( 65 ) is valid and the system of Equation 
( 70 ) is greatly simplified. At small x, j'\/h occurs, since J=2H in 
Solution ( 69 ) and also T-const. One can estimate the terms'with 
derivatives in Equation (70) by introducing the .characteristic wave- 
length X for a disturbance. By comparing the terms with h and dj/dz 
* . • 

in the first of Equation ( 70 ) and with j and d'h/dz in the second of 
Equation ( 70 ), it is easy to see that the terms with derivatives 

j 

always predominate if s . , . 

which is always valid for x<<l. In this case; from the first tv/o 
of Equation (70) with allowance for Equation ( 69 )', it follows that: 
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9 y ’ ' * 2 • (88) 

The constants of integration in Equation (88) w|iich arise during the ' 
solution of Equation (70) have been set equal to zero. They are 
connected with boundary radiation sources which do not depend on the 
lo'cal y displacements. These sources are assumed; to be absent. Sub--, 
stituting Equations ( 69 ) and (88) into the thi'rd of Equation ( 70 ) and 
taking Equation ( 69 ) into account, we shall obtain 

■^’.1 ’ "j' j’?. •• -- 0 , . ■■ 


. - - T ■>. 


(89) 
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Equation (89) is, with allowance for Inequali,ties (87), valid for 
Case (I) at the adiabatic velocity (y*‘ 5/3) and for Case (II) for an 
isothermal v^ (Y®1)- Equation ( 89 ) has an exact ‘solution y a exp(ikz)”. 
Substituting this form of solution into Equation ( 89 ) and introducing * 
the characteristic length from Equation (69), we obtain 


o' - - (kvr^V?. /c*.o. 


Prom this, we have 


i , s? d.j a . it 

The "+" sign is chosen for waves propagating outwards. 


If 


(90) 


(91) 


(92) 


then the gas waves propagate in an isotropic atmosphere, otherwise 
the atmosphere vibrates as a whole (see [19]). Formally the propa- 
gation of some wave follows from Equation ( 86 ) at .w<<v_/2z , but in 

« S 0 

this limiting case the second of Inequalities (. 87 ) 'is violated, and 
therefore. Equation (91) is unsuitable for this case. , If the radia- 
tive damping is sufficiently small 


li ^ 


= I a) 




2) c - Q 

-jf Ky = coT^st Mi 

then, from Equation (91) we have the expansion 
• ^ " 2i. -V* ® c u> J . 


(93) 


(94) 


The first term of Equation (94) determines ‘the increase of amplitude 
in an exponential atmosphere for a constant energy flux. In ful- 
filling the condition Vg<<c3^ the radiative dam’ping in Equation (94) 
is always weaker than the intensification, and .the formation of a 
shock wave occurs. The wave amplitude in the cisse of Equation (94) 
changes, with allowance for Equation ( 69 ), according to the law: 




V. 


c) £al°nl£ti:On of_the_a£ous£l£ £ner£y_f l_ux ^ 


( 95 ) 
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Let U5 find the dependonco of the acoustic enersy flux P on 
optical thickness, using relations that have been derived above. 
Lot us introduce the dimensionless parameters . • 


o^. « -- 5!^ , 


(O' 


O, 


(96) 


We shall adopt from now on the value t»* 2/3 as the- boundary between 
the eransparent and opaque regions. Using Equations (77) and (79) » •* 
we have for the region T>max (2/3, ) • 


F ("c V r ('c,') cx \3 [ 'i • Ct:; 


»/« 


•?/'« 


)l, 


(97) 


'For 2/3 in the region max ( 2/3,Tg) we have from Equation (83) 

(98) 122. 

If To> 2/3 also, then in the region 2/3 <t<t_, using Equation (86), 


we have 


In the transparent region t< 2/3, we have from Equation (95) 

^ “Mj, 

F y) - F (i/^) ( I ^ a; e; C'fc-" 0 . (100) 


s has been shown in Paragraph 5% upon the wave's going over into a 
region with different parameters and for a change^ of the type of wave, 
there is a characteristic quantity which "maintains itself", the 
amplitude of the wave with the smallest damping. In a plane stellar _ 
atmosphere the type of wave changes at the points , Xg and 2/3. ; 
For finding the fraction of acoustic flux which arises at t>>T 2 _ and . 
which flows out into the region x<2/3, we shall use this condition 
of continuity. As follows from the previous examination, three cases 
are possible . , 

I) x^<2/3. ' • ; 


In this case the solution of Equation (*97) immediately goes over 
into Equation (100), the wave velocity and the 'energy flux undergo ah 
abrupt jump of and we obtain ' ‘ : /3^ 

2) ^,>2/3, 110 T^<2/3. ’ - 
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In this case. Solution (97) Joins with. E|quat ion ( 98 ) and, 
moreover, in view of the preservation of the type of wave, the energy 
flux at is continuous, but at t= 2/3, Equation (98) goes over 

into Equation (100) with an abrupt Jump of the acoustic energy flux 
of We have 


F( 2 /a) = ^ [ *5 




(■v;'K.x ■</>), 


( 102 ) 


' 3) Tj>2/3. 

Here a continuous transition from Equation (97) to Equation 
( 98 ) occurs, then a transition with an abrupt Jump of from 
Equation (98) to Equation (99) and,„J!inally , a continuous transition 
from Equation (99) to Equation (100). We have, 

F(Va) -- (f-a 

(103) 

- f e -r )] . 

. d) The out flow of_wav£S_from_a_pla^e_a^mOS2h£re in_the 
presence of a. Corona ' ' ■ 


z 


0 


• Condition (92) is, taking Expressions (68) and ( 69 ) for and 
into account, written down in the form , . . 


60 > 



(IOJ 4 ) 


Here T is the temperature of the gas in a transparent region. As was 
noted in [8], [10] and [11], for turbulent disks ■ accreting onto a 
black hole, the prese.nce of an acoustic energy flux and its transfor- 
mation into heat, and also the heating of the gas by the disk's 
radiation, which creates a non-conservative radiative force, lead to 


the formation of a corona in which the temperature is considerably 
higher than the temperature which is present in an equil- 

ibrium isothermal atmosphere. Coronas also appax^ently exist around . . 
very massive stars. The presence of a corona will lead to the weak- 
ening of Criterion (10^1) for wave transmission, which formally will 

be I’educed to Increases of T and v . This will lead to a broader 

b 

band of outflowing waves and to a stronger acoustical heating, A 

seeming "shining through" of the star's atmosphere occurs upon the /S5 
♦ ■ 

formation of a corona. The increase of the- corona's temperature is 
limited by the fact that the. intrinsic energy losses of the coronal 
gas, which are determined by Bremsstrahlung and inverse Compton radia- ’ 
tion for a radiative flux of photons from the photosphere, become ■ 
significant. A magnetic field can play an Important role in the pro- 
cess of forming a corona (see [12]). . , 

The transition zone from the photosphere from r==2/3 to the corona 
plays an important role In determining the fraction of the acoustic 
energy flux which flows outward. 

If the thickness of the transition zone d„ is much less than , ‘ 

c 

the_ characteristic v/avelength A =2ir/k , k = /v with w from Crlte.r- 

ion (10^1) then one can evidently replace T by T in Criterion (10^). 

c 

In this case the band of frequencies Increases for outflowing waves ; ! 
and the fraction of the outflowing acoustic flux increases. If 

then the presence of a corona must not affect the frequency 
of the outflowing waves and one can use the equilibrium parameters of ■ 
the gas’^’ everywhere. In reality d can be '^JA„; therefore, we shall use’ • 
the equations that have been obtained above with the equilibrium temp-- . 
erature T of the atmosphere for all estimates , 'but ■ let us take into 
account the fact that the frequency of a transmitted wave can be of ' 

the same order or even somewhat less than w. (T) from Inequality (10^)^’ 

C ■ * - 


* With a corona present, waves with oj (T ) <a)<w„ (T) , flow out into the 

corona but their energy flux near tRe base or the corona is reduced , 

by the factor expC^ird /A.,). . . 

c c , • 
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7. Numerical results 


Let us calculate what fraction of the acoustic flux generated 

at large optical depths can flow out into’ the transparent atmosirhere /36 

and heat it. The location of an acoustic wave's transition into a 

shock wave depends on the initial amplitude. We shall consider the 

value of the optical thickness for the transition as the free parameter 

T The calculation has been carried out for two different situations 

out 

in which the role of radiative pressure is great; 1) an accretion 
disk, and 2) the atmosphere of a very massive star. 


In the case of an accretion disk R'x.lOrg, where r^ is the Schwarz- 
schild radius of a black hole, was considered as the region of maxi- •• 
mum energy release. The mass of the black hole' was taken as equal 
to lOM . The turbulence parameter Specifying the parameter '• 

determines all the properties of the disk if a model of the 
vertical structure has been chosen. We considered two variants of the 
vertical structure; the adiabatic structure that has been obtained in ■ 
[8] with a polytropic index n=3 and the n=l structure obtained in [91. 
The results for different frequencies are presented in Table 2. In 
this table the limiting frequency w from Inequality (10^1) is taken 
at the equilibrium atmospheric temperature without allowing for the 
existence of a corona. Therefore, in reality, waves with w slightly . 
lower than can also be transmitted. 

* o * » 


The results of similar calculations for very massive stars 
(see f2]) that are considered as models for quasars are shown in 
Table 3. ' • 

From the results it is evident that the;crude estimate of the 
fraction of the outflowing flux ('v3=P /P ) -that , has been used in [8] • 

g, I • 

is actually valid in order of magnitude for frequencies that are typi- 
cal for a disk (with X'\>z^). The high frequencies rapidly decay ; 
exponentially. fact isolates the characteristic frequency (more'-, 
precisely, the band of frequencies) among the Waves- which heat the- 
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corona of the disk. Let us repeat the factors which lead to such an 
Isolation: 1) an isolated region in the disk' v^lth a maximum temper- ‘ 
ature and maximum energy release;. 2) low frequencies do not p^ss 
through into the corona (besides their generation is also Impeded); 
and 3), the high frequencies decay rapidly. ’ 

8. Comparison with observational data 

* As follows from the previous examination, each object (a disk, 
or a very massive star) possesses a characteristic’ frequency that is. 
determined in Inequality (104). Waves generated by convection with 
a frequency close to w flow out into the atmosphere, disturb it, 
and can lead to the observed fluctuations of brightness. It is of 
interest to compare the characteristic frequencies resulting from 
theory with the observed times for the fluctuations of brightness for 
those objects in which one can suppose accretion disks or very massive 
stars. 

a) The Cy£nus__X-l_X-ray_source 

This source, in which the presence of 'a black hole and accretion 
disk is supposed, is strongly fluctuating in the entire range of 

— '5 

periods from 10 to 10 sec (see [22]). As follows from. Table 2, 
the characteristic fluctuation period caused by. the outflow of 
acoustic waves into the atmosphere amounts to from 5 to 10 millisecs. ' 
The detection of a quasi-period t=;10 milllsec in the x-ray radiation 
of Cygnus X-1 (see [22]) may be associated -with the outflow of waves ■ 
into the, atmosphere above a convective disk. In [22] the existence 
of a quasi-period of vio milllsec is associated with the revolution 
of hot spots around a black hole (see [23]). One can indicate several 
observational differences between these two mechanisms for the fluc- 
tuations. During the revolution of a hot spot .which spirals in toward' 
the black hole, the period of fluctuation in each series must decrease'. 
At the same time the characteristic frequency ‘ which depends on the 
star's mass, does not change as a function of luminosity. If the 
same fluctuations are caused by convection, then' the period in a 


1. <1= 


given series must be approximately constant, but- upon an increase of 
luminosity, the characteristic period increase's approximately pro- 
portionally to the luminosity, as follows from Table 2, Let us 
also notice the difference in .the spectral variation of the two 
mechanisms: the spectrum in a pulse does not change during a revolu- - 

tion, but as a result of convection the hardest spectrum must be at /38 
maximum brightness. However, an increase of the spectrum's hardness ' 
must be noticeable only on the shortest time scales (<10 millisecs). 

Fdr longer intervals the connection of the spectrum with luminosity 
is considerably more complex, since different regions of the disk, 
the flares in which can be uncorrelated, give contributions to : * 

different parts of the spectrum. Therefore, the absence of a simple 
connection of the spectrum with luminosity that has been noted in [24] 
cannot serve as an argument against explaining the spectrum of Cygnus . ; . 
X-1 by a Compton process. Let us notice that the mechanism we have 
considered can give fluctuations that are weakly correlated in time , ^ . 
and can imitate the white noise that has been obtained from an ahalysis; 
of observations of the brightness fluctuations for Cygnus X-1 (see 
[24] and [25]). ‘ • 

b) The nu£lel_of act ive_gal axles and ouasar£- . • 

The nature of quasars and of the compact nuclei of galaxies has . 
been unclear up to now but the physical processes that have been con-',, 
sidered in this paper occur in at least two existing models: the • 

disks around very massive black holes (see, -[3]). and very large stars , 

(see [2]). It is also of interest here to compare the observed pro- ■ • 
perties of variability with the predictions of the model. Rapid flue- . 
tuations of luminosity with a quasi-period of (t,100 days are observed ■ 
in the nucleus of the Seyfert galaxy NGC 4151 ('V'130 days see [26]) 
and in the object OJ 287 ('^184 days), which is -.a' BL Lacertae type 
object (see [27]). These quasi-periods agree well with the character- 
istic [periods of the frequencies in a model for a very massive star 
of '\^10°M^ (see Table 3). According to the data of Table 3, one can .• 
estimate that, in a model of an accretion disk- around a black hole, 
for a hole mass of '^ilO®]yl^ these periods correspond to a luminosity 
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L=0.1 Lq. As is especially noticed by observers (see, for example, 
[27]), "Phase change (which is evidently sharp) with preservatioh 
of the period is a feature of the periodicity of the rapid component 
in the nuclei of Seyfert galaxies". This property agreef.s well with 
the convective-wave hature of rapid fluctuations.’ 

A tendency for an Increase of the quasi-period for the fluctua- 
tions with luminosity, an approximate constancy • of the period in 
each given series of observations, and also an Increase of the hard- 
ness of the spectrum at maximum brightness are the general properties 
of the convective-wave mechanism for brightness as fluctuations both 
in a model for a very large ma.ss star., and also in a model for a turbu 
lent convective disk around a very massive ’black hole. 

However, the role of radiation pressure Is' very great and the 
damping of acoustic waves is very strong in both these models.. It is 
clear from Table 3 that in this case the outflowing acoustical flux 
amounts to no more than '\> 1 % of the flux generated at a large optical 
depth, which is significantly less than the observed amplitude of 
variability. It appears to us that the transfer of energy by other 
types of waves (for example, by magneto-acoustic waves) will also be 
strongly impeded in this case, since damping of them by radiative 
friction will occur effectively. This difficulty becomes especially 
important in a model of a very massive star, where there the mechan- 
isms of variability that are specific for an accretion disk are not 
present (also see our paper [28]). 

In conclusion, let us quote the main results of the paper. 

1. Equations have been obtained which take into account v/c 
terms and which describe the propagation of plane linear waves in a 
medium with high radiative pressure, wl^iere 'the • radiative transfer is 
described by the Eddington approximation. . • 


32 




2. The damping and transformation of wayes have been invest!- 

I 

gated in model problems (an abrupt Jump of opacity and a zone of 
variable opacity). • ■ 

3. The outflow is considered of waves’ into the atmosphere and 
corona of an accretion disk and of a very massive star.' , 

4. The characteristic times of variability have been obtained 
for the Cygnus X-1 source and for the active nuclei of galaxies and 
of quasars and the difficulties of the model for a very massive star* 
in the latter case have been indicated. 

The authors are grateful to Ya. B. Zel'dovich for important 
remarks. We thank Ch. Kunash and M. M. Basko for making available 
the material of [21]. 




TABLE 1 , 

* t 

The complex amplitudes of waves which were formed in the transition 
region with variable opacity. The amplitude of an incident radiation ’ . 
wave is A^^l, the increase of a photon path Cp^/C^ s 10^, and the velo- , 
city ratio of the gas wave to the radiation wave is v /v =0.1. The 

,s ^ * 

frequency of the wave w is selected such that 3 , w/o)^=0,l, ' • 

and w/wg=10. The parameter ^ is the half-width of the transition 
zone (see Equation (62)), A^ and Ag are the amplitudes of the reflected 
waves at the point , and and Bg are those of the transmitted 

waves at the point ^ = . The case^/Aj = Cis an analytical calcula- 

tion of the abrupt jump (see Equation (59)) » and the remaining variant’s, 
are by numerical calculation. 



A, 

A 3, 

6. 

' 


0 

-0.247 + 0.247 i 

0.553 -’0.447 i 

I.X05 0 i 

0.045 + 0.155 i 1 

O.OI 

-0.248 + 0.256 i. 

0.478 - 0.2,00 i 

J..I79 - 0.09C- ^ 

0.052, + 0.139 ( 

O.I 

-0.250 t 0.229 i 

0.47'8 - 0.234 ( 

1. 169 + 0.091 'i 

0.051 .p 0.139 i 

1 I.O 

-0.165 + 0.026 i 

0.494 - 0.2a3 ^ 

-0.027 + X.022 \ 

0.040 + 0.I3S V 

1 10 

1 

-0.0052- 0.0012 V 

0«3X6 + 0»XVx 

-0.056 -p 0.022 ^ 

-0.049 + 0.0S3V 



5 
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3^1 
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The outflowing acoustical flux P in an accretion disk around a 
black hole of M=10M^ for an outflow optical thickness (meaning; 

outflow optical thickness) for different values of the parameters m= 

= ^g/^r» wavelengths X relative to the half-thickness of the . ‘ 
disk 0 Q with corresponding values of the period t, and the oscilla- 
.tion frequencies w in fractions of the limiting frequency w obtained • 
without allowance for the heating of the atmosphere. 


rn 


> 

u) 

•t 

•r: 


r/F. 


j 

sec . 1 

■t. "A5 
- out 

X.. ' qu 

‘A rO.l 
-out „ 

O.ObI 

0.010 

O'. I 


O.I 

1 

i 

2.2 

4.4 

5.0 

2.5 

II 

22 

1.6 

3.3 

0.007 

0.019 

0.006 

0.019 

[M 

0.27 

i 

a 

i. 

<> 

0.0 

2.4 

5.1 

1.7 

II 

33 

4.5 

13 

0.23 

0.045 

0.22 

0.044 

0.20 

0.043 

0.3 

b] 

o.oi 

I 

1 

X 

4.2 

0.4 

JO 

5 

3.7 

i.9 

0.06 

0.03 

0.036 

0.023 

0.029 

0.010 

0.017 

0.010 

(el 

0.032 

1 

X 

2.7 

5 

3.8 

0.18 

0.040 

0.037 

0.030 
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The outflowing acoustical flux in a very massive star as a 
quasar model for different radii R (in units of the Schwarzschlld 
radius r^) , The notations are the same as in Table 2. 
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8,5-10“® 

100 


10 

8.5-10“'* 

.100 


10 


I.O'IO'^ 4.7 O.OS 
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